Quantum mechanical time-delay matrix in chaotic scattering 
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We calculate the probability distribution of the matrix Q — —ihS~ dS/dE for a chaotic system 
with scattering matrix S at energy E. The eigenvalues Tj of Q are the so-called proper delay times, 
introduced by E. P. Wigner and F. T. Smith to describe the time-dependence of a scattering process. 
The distribution of the inverse delay times turns out to be given by the Laguerre ensemble from 
random-matrix theory. 
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Eisenbud and Wi gner IS] introduced the notion of 
time delay in a quantum mechanical scattering problem. 
Wigner's one-dimensional analysis was generalized to an 
N x N scattering matrix S by Smith ||, who studied 
the Hermitian energy derivative Q — — ifiS~ 1 dS / dE and 
interpreted its diagonal elements as the delay time for a 
wave packet incident in one of the N scattering channels. 
The matrix Q is called the Wigner-Smith time-delay ma- 
trix and its eigenvalues n , t%, . . . , tjv are called proper 
delay times. 

Recently, interest in the time-delay problem was re- 
vived in the context of chaotic scattering p|. There is 
considerable theoretical and experimental |s|-|l0t ev- 
idence that an ensemble of chaotic billiards containing a 
small opening (through which N modes can propagate at 
energy E) has a uniform distribution of S in the group 
of N x N unitary matrices — restricted only by funda- 
mental symmetries. This universal distribution is the cir- 
cular ensemble of random-matrix theory [11], introduced 
by Dyson for its mathematical simplicity [12]. The eigen- 
values e l< ^ n of S in the circular ensemble are distributed 
according to 
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with the Dyson index /? = 1, 2, 4 depending on the pres- 
ence or absence of time-reversal and spin-rotation sym- 
metry. 

No formula of such generality is known for the time- 
delay matrix, although many authors have worked on this 
problem P,|l3|-p3|| . An early result, (trQ) = th, is due 
to Lyuboshits |13[ , who equated the ensemble average of 
the sum of the delay times tr Q = J^n=i T « to the Heisen- 
berg time tjj — 2irh/A (with A the mean level spacing 
of the closed system) . The second moment of tr Q was 
computed by Lehmann et al. p8| and by Fyodorov and 
Sommcrs |l9]. The distribution of Q itself is not known, 
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The trace of Q determines 
and is therefore sufficient for 
For applications 



the density of states |2 
most thermodynamic applications |21^ 
to quantum transport, however, the distribution of all 
individual eigenvalues t„ of Q is needed, as well as the 
distribution of the eigenvectors (2f| . 

The solution of this 40 year old problem is presented 
here. We have found that the eigenvalues of Q are inde- 
pendent of S J2(| . The distribution of the inverse delay 
times 7„ = l/r„ turns out to be the Laguerre ensemble 
of random- matrix theory, 
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but with an unusual iV-dependent exponent. (The func- 
tion P is zero if any one of the t„'s is negative.) The 
correlation functions of the r„'s consist of series over 
(generalized) Laguerre polynomials p7[ , hence the name 
"Laguerre ensemble" . The eigenvectors of Q are not in- 
dependent of S, unless f3 = 2 (which is the case of broken 
time-reversal symmetry). However, for any j3 the corre- 
lations can be transformed away if we replace Q by the 
symmetrized matrix 
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which has the same eigenvalues as Q. The ma- 
trix of eigenvectors U which diagonalizes Qe — 
f/diag(ri, . . . ,tn)W is independent of S and the t„'s, 
and uniformly distributed in the orthogonal, unitary, or 
symplectic group (for = 1, 2, or 4, respectively). The 
distribution @ confirms the conjecture by Fyodorov and 
Sommers |l^] that the distribution of tr Q has an alge- 
braic tail oc (trQ)~ 2 -W 2 . 

Although the time-delay matrix was interpreted by 
Smith as a representation of the "time operator" , this 
interpretation is ambiguous Jh|. The ambiguity arises 
because a wavepacket has no well-defined energy. There 
is no ambiguity in the application of Q to transport 
problems where the incoming wave can be regarded 
monochromatic, like the low-frequency response of a 
chaotic cavity |2l[|22|,|28|| or the Fermi-energy dependence 
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of the conductance p5[ . In the first problem, time de- 
lay is described by complex reflection (or transmission) 
coefficients R mn (ui), 



Rmn{u) — R 
-Rmn(O) = S. 
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lUJT rnn + 0(u J 2 )}, (4a) 
= ImhS m 1 n dS mn /dE. (4b) 



The delay time r mn determines the phase shift of the a.c. 
signal and goes back to Eisenbud [Q. With respect to 
a suitably chosen basis, we may require that both the 
matrices R m n(fi) and r mn are diagonal. Then we have 



-RmnM = 6 mn [l + iujT m + 0{u> 2 )}, 



(5) 



where the r m (m = 1, . . . , N) are the proper delay times 
(eigenvalues of the Wigner-Smith time-delay matrix Q). 
For electronic systems, the O{to) term of R mn {u>) is the 
capacitance. Hence, in this context, the proper delay 
times have the physical interpretation of "capacitance 
eigenvalues" f2S[] , 

We now describe the derivation of our results. We start 
with some general considerations about the invariance 
properties of the ensemble of energy-dependent scatter- 
ing matrices S(E), following Wigner PQ] , and Gopar, 
Mello, and Biittiker The N x N matrix S is uni- 

tary for (3 = 2 (broken time-reversal symmetry), uni- 
tary symmetric for (3=1 (unbroken time-reversal and 
spin-rotation symmetry), and unitary self-dual for /3 = 4 
(unbroken time-reversal and broken spin-rotation sym- 
metry). The distribution functional P[S(E)] of a chaotic 
system is assumed to be invariant under a transformation 



S(E) -► VS{E)V, 
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where V and V' are arbitrary unitary matrices which do 
not depend on E (V = V T for (3 = 1, V = V R for 
(3 = 4, where T denotes the transpose and R the dual 
of a matrix) . This invariance property is manifest in the 
random-matrix model for the ^-dependence of the scat- 
tering matrix given in Ref. p2|| . A microscopic justifica- 
tion starting from the Hamiltonian approach to chaotic 
scattering ]3l| is given in Ref. |32|. Eq. (||) implies with 
V = V = iS- 1 ' 2 that 



P(S, Qe) = P(-l, Qe) 



(7) 



Here P(S, Qe) is the joint distribution of S and Qe, de- 
fined with respect to the standard (flat) measure cLQe for 
the Hcrmitian matrix Qe and the invariant measure dS 
for the unitary matrix S. From Eq. (Q) we conclude that 
S and Qe are statistically uncorrelated; Their distribu- 
tion is completely determined by its form at the special 
point S = — 1. 

The distribution of S and Qe at S = — 1 is com- 
puted using established methods of random-matrix the- 
ory Jll]|n| ■ The N x N scattering matrix S is expressed 
in terms of the eigenvalues E a and the eigenfunctions 



Vw of the M x M Hamiltonian matrix H of the closed 
chaotic cavity H, 
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The Hermitian matrix Ti. is taken from the Gaussian or- 
thogonal (unitary, symplectic) ensemble ]Tl"| ], P(H) oc 
exp(-/37r 2 trW 2 /4A 2 M). This implies that the eigenvec- 
tor elements ij)j a are Gaussian distributed real (complex, 
quaternion) numbers for (3 = 1 (2, 4), with zero mean and 
with variance M~ x , and that the eigenvalues E a have 
distribution 



P{{E a })^\{\E,-E v \ P X{< 



-f3TT 2 El/4A 2 M 
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The limit M — > oo is taken at the end of the calculation. 

The probability P(— 1, Qe) is found by inspection of 
Eq. (|sj) near S = — 1. The case S 1 = — 1 is special, be- 
cause S equals —1 only if the energy E is an (at least) 
./V-fold degenerate eigenvalue of H. For matrices S in a 
small neighborhood of —1, we may restrict the summa- 
tion in Eq. (||) to those TV energy levels E a , a = 1, . . . , N, 
that are (almost) degenerate with E (i.e. \E — E a \ <C A). 
The remaining M — N eigenvalues of H do not contribute 
to the scattering matrix. This enormous reduction of the 
number of energy levels involved provides the simplifica- 
tion that allows us to compute the complete distribution 
of the matrix Qe- 

We arrange the eigenvector elements ip na into anNxN 
matrix = ^M 1 / 2 . Its distribution P(^) oc 

exp(— (3tr /2) is invariant under a transformation 
^ — > "JO, where O is an orthogonal (unitary, symplec- 
tic) matrix. We use this freedom to replace $ by the 
product 'SO, and choose a uniform distribution for O. 
We finally define the N x N Hermitian matrix Hij = 
J2a=i Oia{E a — E)Oj a . Since the distribution of the en- 
ergy levels E a close to E is given by ]~L<i/ I-^m — 
[cf. Eq. (g)], it follows that the matrix H has a uniform 
distribution near H = 0. We then find 



5 = -1 + (iTB/h)^- 1 ^- 1 , 
Q e =t h ^-H-\ 
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Hence the joint distribution of S and Q^i at S = — 1 is 
given by 

P(-1,Q B )« / d*d J ffe- /3tr *' I ' t / 2 

x (5(* t - 1 i/*- 1 ) (5(Q £ - r H *t-i^-i) 

xSiQE-TH^- 1 *- 1 ). (11) 

The remaining integral depends entirely on the positive- 
definite Hermitian matrix T = ty^. In Refs. pTJ and 
[p3l it is shown that 
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J = /" dr (det r) (/3 ~ 2)/2 /(r)e(r), (12) 

where 0(r) = 1 if all eigenvalues of T are positive and 
otherwise, and / is an arbitrary function of T = \I/^E 
Integration of Eq. @ with the help of Eq. @ finally 
yields the distribution (Q) for the inverse delay times and 
the uniform distribution of the eigenvectors, as adver- 
tised. 

In addition to the energy derivative of the scattering 
matrix, one may also consider the derivative with respect 
to an external parameter X, such as the shape of the 
system, or the magnetic field Jl^j2^]. In random- matrix 
theory, the parameter dependence of energy levels and 
wavefunctions is described through a parameter depen- 
dent M x M Hermitian matrix ensemble, 



H(X) = H + M- 1/2 XH', 



(13) 



where TL and TL' are taken from the same Gaussian 
ensemble. We characterize dS/dX through the sym- 
metrized derivative 



Q 



x 



jO-l/2 dS q-1/2 

%b dx ' 
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by analogy with the symmetrized time-delay matrix Q e 
in Eq. (|J). To calculate the distribution of Qx, we 
assume that the invariance (0) also holds for the in- 
dependent ensemble of scattering matrices. (A random- 
matrix model with this invariance property is given in 
Ref. p4|). Then it is sufficient to consider the special 
point S = -l. From Eqs. fllObl ) and |[|) we find 

Qx = * t " 1 -ff'* _1 , P(H') oc exp(-/3tr/f' 2 /16), (15) 



where H' 



r ip %ip jv - Acalcu- 



-(r H /h)M-y^ 
lation similar to that of the distribution of the time-delay 
matrix shows that the distribution of Qx is a Gaussian, 
with a width set by Qe, 



P(S,QE,Qx)^(detQ E r 2 ^ 3+;i ^ 2 
x exp 



f tr (thQe 1 + 1 -(t h Q- e x Q x ) 2 



(16) 



The fact that delay times set the scale for the sensitiv- 
ity to an external perturbation in an open system is well 
understood in terms of classical trajectories | p5[ , in the 
semiclassical limit N — > oo. Eq. (16) makes this pre- 
cise in the fully quantum-mechanical regime of a finite 
number of channels TV. Correlations between parameter 
dependence and delay time were also obtained in Refs. 
iP^^ot , for the phase shift derivatives d(j)j/dX. 

In summary, we have calculated the distribution of 
the Wigner-Smith time-delay matrix for chaotic scatter- 
ing. This is relevant for experiments on frequency and 
parameter-dependent transmission through chaotic mi- 
crowave cavities J9|,[l0| or semiconductor quantum dots 



with ballistic point contacts |3(J. The distribution ([!]) 
has been known since Dyson's 1962 paper as the circular 
ensemble jl2| . It is remarkable that the Laguerre ensem- 
ble (0) for the (inverse) delay times was not discovered 
earlier. 
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